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ABSTRACT 


This  report  summarizes  the  activities  carried  out  under  Signal 
Corps  Contract  No.  DAAB07-68-C-0381  (Texas  A  &  M  Research  Foundation 
Project  586)  during  the  contract  period  15  December  1968  through 
15  June  1969. 

During  this  period  a  set  of  equations  which  incorporate  the 
momentum  and  energy  sources  and  sinks  due  to  the  trees  and  foliage 
were  developed  for  analog  simulation  of  the  forest  atmosphere. 

The  region  above  and  within  the  forest  are  treated  separately  but 
are  coupled  at  the  top  of  the  canopy. 
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I.  Introduction 

The  trees  and  attendant  foliage  of  the  forest  provide  a  complex 
distribution  of  sources  and  sinks  of  heat,  moisture  and  momentum. 

In  view  of  the  diverse  characters  of  the  atmospheric  layers  above 
and  interior  to  the  forest,  the  forested  boundary  layer  is  most 
conveniently  considered  in  two  sections:  the  upper  free-air  region 
and  the  lower  forest  region. 

Experimental  systems  of  equations,  with  coupling  effected  through 
energy  and  momentum  considerations  at  the  top  of  the  forest  canopy, 
have  been  developed  to  simulate  these  two  regions.  The  systems  are 
broadly  based  on  the  LLMM  with  some  significant  modifications. 

These  sets  of  equations  are  now  being  programmed  for  analog  solution. 

The  equations  used  in  the  simulation  are  listed  in  the  Appendix 
anu  are  discussed  in  Section  II.  Some  of  the  details  of  the  preliminary 
simulation  plan  are  given  in  Section  111. 
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II.  Defining  Relationships 

The  equations  developed  to  describe  the  forested  boundary  layer 

are  based  on  the  assumptions  listed  in  Semi-Annual  Report  1  (hereafter 

referred  to  as  SRI)  plus  additional  specifications  for  the  free  air- 

canopy  interface  and  the  vertical  distribution  of  the  exchange 

coefficient  for  momentum.  K  .  The  conditions  at  the  interface  are 

m 

characterised  by  the  conservation  of  energy  at  the  top  of  the  canopy. 

The  momentum  exchange  coefficient  is  assumed  constant  in  the  ground 

layer  (depth  A')  and  in  the  air-canopy  layer,  (d  -  A)  4  z  $  (d  +  A) , 

where  d  is  the  height  of  the  forest.  The  values  of  K  for  these  layers 

m 

are  based  on  Deacon's  Beta-equation  calculated  independently  from  the 
vertical  gradients  of  temperature  and  wind  speed  just  above  the  forest 
floor  and  just  above  the  top  of  the  canopy. 

A.  Free  Air  Section 

Equations  (1)  to  (18)  in  the  Appendix  define  the  temporal 

-  riations  of  momentum,  temperature  and  vapor  pressure  in  the  free 

air  above  the  forest.  These  are  essentially  the  same  as  the  atmospheric 

section  of  the  LLMM  with  the  exception  of  equations  (14)  through  (16) . 

In  view  of  the  current  uncertainties  in  the  magnitude  of  the  eddy 

diffusion  coefficients,  they  are  assumed  equal  in  equation  (14). 

The  vertical  distribution  of  the  momentum  exchange  coefficient 

above  the  air-canopy  layer  (Equation  (15))  is  the  same  as  that  in  the 

LLMM  except  that  the  entire  distribution  is  lifted  a  distance  equal 

to  the  height  of  the  forest,  d.  The  exchange  coefficient,  K  ,  govs  to 

m 
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zero  «t  the  top  of  the  boundary  layer  plus  the  increment  d  rather  than 

at  Z,  the  top  of  the  boundary  layer  Itself.  Thus  the  friction 

layer  is  deeper  than  it  would  be  over  a  smooth  surface  by  an  amount 

equal  to  the  height  of  the  forest.  This  modification  was  selected 

from  one  of  several  possibilities  because  it  retained  the  original 

vertical  distribution  yet  took  into  consideration  the  fact  that  the 

roughness  was  greatly  increased.  The  coefficient  b,  in  Equation  (15) 

is  determined  in  Equation  (16)  from  the  requirement  that  K  be  continuous 

in 

at  the  boundary  between  the  air-canopy  layer  and  the  free  air  region. 

It  is  recognized  that  this  approach  does  not  eliminate  potential 

difficulties  with  poles  in  the  K  ..  .  . ,  solution  as  well  as  some 

m(d  +  A) 

l*;ck  of  reality,  the  latter  because  the  Deacon  3  equation  was  developed 
from  data  collected  in  non-forested  regimes.  However,  an  alternative 
is  under  consideration,  if  these  potential  difficulties  become  real. 

B.  Free  Air-Canopy  Section 

The  relationships  for  the  air-canopy  layer  are  given  in  Equations 
(19)  through  (38)  of  the  Appendix.  The  most  basic  equations  are  (19) 
and  the  combination  (32),  (33)  which  define  respectively,  the  energy 
balance  and  eddy  stresses  at  canopy  top. 

The  energy  balance  (Equation  (19))  recognized  that  heat  may  be 
transferred  away  from  the  interface  by  eddy  fluxes  of  heat  and  moisture 
downward  through  the  canopy  as  well  as  upward  through  the  free  air.  The 
contribution  from  plant  metabolism  is  much  smaller  than  the  other 
factors  and  has  been  neglected. 
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The  net  radiation  is  calculated  as  in  the  LLMM  except  for  the 
modification  of  the  short  wave  transmissivity  of  the  atmosphere 
mentioned  in  SRI.  In  addition  the  reflection  of  long  wave  radiation 
by  the  canopy  is  neglected,  since  it  represents  less  than  2  or  3 
percent  of  the  up-going  long  wave  radiation. 

The  eddy  fluxes  (Equations  (25)  through  (28))  are  standard  in 
form.  The  parameter,  £,  in  the  equation  for  the  vapor  pressure  at 
level  d  (Equation  (29))  is  related  to  the  conductivity  of  the  canopy 
which  is  a  function  of  a  number  of  variables  including  leaf  area  and 
stomatal  resistance.  It  thus  is  a  problem  variable  dependent  on 
canopy  type  and  season.  Studies  currently  underway  are  seeking 
methods  by  which  this  variable  may  be  parameterized. 

The  stresses  at  the  canopy  top  are  calculated  on  the  implicit 
assumption  of  nearly  linear  change  of  wind  with  height  through  the  layer 
d-A<z<d+A  (Equation  (32),  (33)).  This  necessarily  puts 
restrictions  on  the  depth  A,  which  will  depend  on  the  height  and 
type  of  forest.  The  momentum  exchange  coefficient  has  been  assumed 
constant  through  this  layer ,  and  is  specified  by  the  value  at  d  +  A, 
(from  Equation  (37))  on  the  assumption  that  the  Deacon  profile  describes 
the  wind  in  the  layer  immediately  above  the  canopy  top.  The  roughness 
parameter,  zq,  for  the  forest  is  determined  using  the  empirical 
relationship  developed  by  Kung  (1961).  Recent  wind  tunnel  studies 
of  model  foreBts  (Hsi  and  Nath,  1968)  indicate  that  this  relationship 
may  be  appropriate. 
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C.  The  Forest  Section 


The  forest  region  is  naturally  the  most  difficult  to  formulate 

not  only  because  of  the  complex  vertical  distribution  of  sources 

but  also  because  the  source  terse  themselves  are  complicated  and 

incompletely  understood.  The  momentum  sink  of  the  forest  is  Included 

in  the  momentum  equations  (If)  and  (2f)  by  the  form  drag  termB 

F  and  F  defined  In  Equations  (9f)  and  (lOf).  These  are  assumed 
x  y 

constant  throughout  the  depth  of  the  forest  but  can  be  made  a  function 
of  height  if  it  appears  necessary  at  a  later  time  through  the 
representative  coverage  factor  A.  The  coverage  factor  A  is  a  characteristic 
of  the  forest  and  is  the  effective  frontal  area  per  unit  depth  of  forest 
per  unit  horizontal  area.  The  drag  coefficient,,  Cp,  is  the  local  drag 
coefficient  for  the  entire  forest  depth. 

The  temporal  change  of  temperature,  Equation  (3f ) ,  can  be  viewed 
as  the  sum  of  a  number  of  contributions:  horizontal  advection, 
convective  flux  divergence,  radiation  flux  divergence,  and  plant 
metabolism.  The  last  factor  is  small  relative  to  the  others  and  can 
be  neglected  as  it  is  in  the  energy  balance  at  the  canopy  top.  The 
sensible  heat  (convective)  fluxes  are  calculated  in  the  usual  way 
(Equation  (Ilf),  (12f)).  There  is  some  difficulty  however  in  determining 
the  proper  relationship  for  R,  the  radiation  flux  divergence.  This 
factor  must  now  include  the  radiation  from  the  trees  and  foliage  as 
well  as  the  radiation  flux  from  the  atmospheric  water  vapor.  There 
is  some  evidence,  e.g.  Denmead,  1964,  that  the  source-sink  distributions 
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of  radiant  energy  are  not  uniform,  so  that  a  simple  constancy  of 
radiation  flux  divergence  through  the  depth  of  the  forest  may  not  be 
adequate.  Studies  are  currently  under  way  to  develop  a  computational 
definition  of  the  radiation  term  based  on  atmospheric  variables. 

The  moisture  expression,  Equation  (4f)  includes  the  contribution 

to  the  moisture  content  of  the  air  made  by  evapo-transpiration.  A 

defining  expression  for  this  term,  M,  is  being  sought  in  current  studies 

of  forest  energy  budgets.  The  stress  and  the  convective  and  evaporative 

flux  equations  (7f),  (8f) ,  (Ilf),  and  (12f ) ,  are  standard.  The  eddy 

diffusion  coefficients  are  assumed  equal  as  they  are  above  the  forest 

(Fcjation  (16f))  but  the  defining  expression  for  is  different.  The 

defining  equation  for  K  given  in  Equation  (17f)  follows  from  the 

assumption  of  a  linear  distribution  between  the  ground  section  and  the 

9h 

ai  -canopy  layer.  The  pressure  gradient,  ,  is  considered  constant 

through  the  depth  of  the  forest  and  equal  to  the  surface  pressure  gradient. 

The  assumption  of  zero  vertical  velocity  is  tentatively  maintained 
pending  first  testing  results  and/or  more  definitive  wind  data  for 
forested  regions.  This  assumption  also  holds  for  the  Forest  Surface 
Section  discussed  below. 

D.  The  Forest  Surface  Section 

The  conditions  in  the  forest  surface  layer,  from  one  meter  below 
the  ground  surface  to  a  height  A'  above  the  surface  are  specified  by 
conservation  of  energy  at  the  surface  and  the  assumption  of  constancy 
of  the  momentum  flux  through  the  layer.  The  energy  budget  at  the 
forest  surface  given  in  Equation  (18f)  is  the  same  as  that  used  in 
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Che  LLMM.  The  net  radiation  considers  the  down-coming  and  up-going 
radiation  in  both  the  short  wave  and  long  wave  parts  of  the  spectrum. 

The  down-coming  short  wave  radiation  is  expressed  as  some  fraction  of 
that  arriving  at  the  top  of  the  canopy,  that  fraction  given  by  the 
transmissivity  (Equation  <20f)).  As  in  the  air-canopy  section,  the 
reflected  long  wave  radiation  is  neglected  (Equation  (23f)).  The  other 
expressions  for  the  forest  surface  layer  (Equation  (24f)  through  (37f)) 
have  already  been  discussed  in  SRI,  except  that  the  depth  of  the  surface 
layer  has  been  expressed  as  the  variable  A*  rather  than  as  a  constant  8m. 
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III.  Layers  and  Levels 

The  choice  of  layers  and  levels  of  Che  forest  simulator  Is 
somewhat  arbitrary  for  the  free-air  section.  Within  the  forest,  however, 
these  must  depend  on  the  depth  of  the  forest.  Initially  the  forested 
boundary  layer  will  be  treated  approximately  as  follows  (where 
d  is  the  height  of  the  forest  and  A'  and  A  define  the  forest  surface 
section  and  the  air-canopy  section,  resp.) 

1,  Free  Air  Section  (Equations  (1)  through  (18)) 


Levels  (m) 

Layers  (m) 

(a) 

x  m  d  +  A 

d  to  d  +  26 

(b) 

x  -  d  +  (32  -  A) 

d  +  2A  to  d  +  64 

(O 

x  -  d  +  107 

d  +  64  to  d  +  150 

(d) 

x  «  d  +  200 

d  +  150  to  d  +  250 

Free  Air- 

-Canopy  Section  (Equations  (19)  through  (38)) 

Levels 

Layers 

z  ■  d 

d  -  A  to  d  +  A 

Forest  Section 

Levels 

Layers 

(a) 

A 

surface  to  2A' 

(b) 

§  -  (A  -  A') 

2A’  to  d  -  2A 

(c) 

d  -  A 

d  -  24  to  d 

4.  Forest  Surface  Section 

Levels  Layers 

surface  -lm  to  A' 
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IV.  Experimental  Simulation 

The  equations  are  being  programed  for  solution  in  the  analog 
computer.  In  this  initial  simulation  the  top  of  the  canopy  is 
assumed  to  be  40a  and  4  and  4'  are  5  and  l.Sa  resp. 

The  layers  and  levels  there,  are  as  follows. 


Layer 

Level 

300  -  200m 

250m 

200  -  100 

150 

M 

a 

100  -  50 

75 

V 

u 

U* 

Lr 

O 

o 

45 

40  m 


4J 

40  -  30 

35 

m 

0) 

£ 

o 

30-3 

15 

Pn 

o 

1 

1. 

One  of  the  primary  objectives  of  the  initial  simulation  will  be 
the  test  of  the  energy  and  momentum  source  terms  proposed  for  the 
forest  section.  In  addition,  studies  will  be  made  to  determine  how 
sensitive  the  model  is  to  other  variables  and  problem  parameters. 
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GLOSSARY 


A 

characteristic  tree  coverage 

(cm  '*) 

a 

ratio  of  molecular  weight  of  water  to 
molecular  weight  of  dry  air 

(non-dimensional) 

CD 

local  drag  coefficient 

(non-dimensional) 

C 

P 

specific  heat  of  air  at  constant  pressure 

(cal/gm  deg) 

d 

height  of  the  effective  top  of  the  canopy 

(cm) 

V 

Integral  exchange  coefficient  for  momentum 
between  surface  and  height  A’ 

(cm/Bec) 

e 

mean  vapor  pressure 

(mb) 

ef 

representative  vapor  pressure 

(mb) 

*0,8 

surface  saturation  vapor  pressure 

(mb) 

ed,s 

saturation  vapor  pressure  at  the  temperature  of 
the  canopy  top 

(mb) 

F 

c 

cloud  factor  for  insolation 

(non-dimensional) 

f 

Coriolis  parameter 

(rad/sec) 

F  ,F 
x*  y 

components  of  drag  force  due  to  trees  and 
foliage 

(cm/sec) 

G 

thermal  resistance  of  surface  litter 

2 

(cm  sec  deg/cal) 

g 

acceleration  due  to  gravity 

(cm/sec2) 

H 

hour  angle  (zero  for  local  apparent  noon) 

(rad) 

h 

height  of  a  constant  pressure  surface 

(cm) 

1 

mean  solar  constant 

2 

(cal/cm  sec) 

1 

index 

j 

albedo 

(noo-dirnensiorial) 
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. d 


*4. 

exchange  coefficient  for  heat 

2 

(cm  /sec) 

K 

B 

exchange  coefficient  for  momentum 

(cm^ /sec) 

K 

▼ 

exchange  coefficient  for  water  vapor 

(ca?  /sec) 

k 

.40  (Von  Kansan’ a  constant) 

(non -dimens ional) 

L 

latent  heat  of  vaporisation  of  water 

(cal/gm) 

h 

net  longwave  radiation 

2 

(cal/cm  sec) 

m 

enplrical  radiation  factor 

(non-dimensional) 

n 

empirical  radiation  factor 

(riT1/2) 

M 

moisture  source  tern  for  forest 

(mb/sec) 

N 

turbidity 

(non -dimens ional) 

P 

86,400  (diurnal  period) 

(sec) 

P 

atmospheric  pressure 

(mb) 

Q 

energy  addition  per  unit  mass  from 
non-adiabatic  processes 

(cal/gm) 

q 

specific  humidity 

(non-dimensional) 

qc 

convective  heat  flux,  positive  upward 

2 

(cal/cm  sec) 

qc,d+ 

convective  heat  flux  in  air  layer  just  above 
canopy,  positive  upward 

2 

(cal/cm  sec) 

"c.d- 

convective  heat  flux  in  canopy  just  below  the 
level  d,  positive  downward 

2 

(cal/cm  6ec) 

qc 

evaporative  heat  flux,  positive  upward 

2 

(cal/cm  sec) 

qe,d+ 

evaporative  heat  flux  in  air  layer  just  above 
canopy,  posicive  upward 

2 

(cal/cm  sec) 

Vd- 

evaporative  heat  flux  in  canopy  just  below 
the  level  d,  positive  downward 

2 

(cal/cm  sec) 
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R 

R 

a 

*N 

S 


s 

o 

T 


T' 


z 

z 

o 

zero 
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6 

A 


soil  heat  flux,  positive  downward 
radiational  cooling  or  wanning 
gas  content  for  dry  air 
net  radiation 
wind  speed 

net  shortwave  radiation 

surface  moistness 

mean  temperature  of  air 

mean  soil  temperature 

representative  soil  temperature 

time 

mean  east-west  component  of  wind 

mean  north-south  component  of  wind 

geos trophic  wind  components 

east-west  coordinate,  positive  eastward 

north-south  coordinate,  positive  northward 

105,000 

vertical  coordinate,  positive  upward 

surface  roughness  length 

as  a  subscript  indicates  value  at  z  *  0, 
excepting  sq  and 

stability  parameter  in  the  Deacon  profile 

solar  declination 

1/2  depth  of  air-canopy  layer 


(cal/cm  sec) 

(deg/sec) 

2  2 

(cm  /sec  deg) 

2 

(cal/cm  sec) 

(cm/  sec) 

2 

(cal/cm  sec) 

2 

(cal/cm  sec  mb) 
(deg  C.) 

(deg  C.) 

(deg.  C.) 

(sec) 

( cm/sec) 
(cm/sec) 
(cm/sec) 

(cm) 

(cm) 

(cm) 

(cm) 

(cm) 


(non-dimensional) 

(rad) 

(cm) 
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A* 

depth  of  forest  surface  layer 

(cm) 

Ah 

vertical  interval  between  adjacent  simulator 
levels 

(cm) 

£ 

emissivity 

(non-dimensional) 

solar  zenith  angle 

(rad) 

e 

?(.'an  potential  temperature 

(deg  C.) 

X 

volumetric  heat  capacity  of  soil 

3 

(cal /cm  deg) 

V 

thermal  conductivity  of  soil 

(cal/cm  sec  deg) 

e 

moisture  parameter  for  the  canopy 

(cal/cm  sec  nib) 

i 

3 • 14 *  .  •  • 

(non-dimensional) 

p 

air  density 

(gm/cm3) 

0 

5.67  x  10~5  Stefan-Boltzmann  constant 

O  a| 

(ergs/cm  sec  deg  ) 

T 

x 

component  of  *  in  x  direction 

2 

(dynes/ cm  ) 

T 

v 

component  of  *  in  /  direction 

2 

(dynes /cm  ) 

J 

* 

latitude 

(deg) 

X 

forest  transmissivity 

(non-dimensional) 

1» 

solar  distance  factor 

(non -dimens  tonal) 

bi 

7.3  x  10”5  (angular  velocity  of  earth's 
rotation) 

(rad/sec) 
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